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Abstract 

The formulas for the determination of the coefficients of the spherical harmonic expansion of 
the disturbing potential of the earth are defined for data given on a sphere. In order to determine 
the sDhericaf harmonic coefficients, the gravity anomalies have to be analytically downward 
continued from the earth's surface to a sphere - at least to the ellipsoid. The goal of tins paper is to 
continue the gravity anomalies from the earth’s surface downward to the ellipsoid using recent 
elevation models. The basic method for the downward continuation is the gradient solution (the gi 
term). The terrain correction has also been computed because of the role it can play as a correction 
term when calculating harmonic coefficients from surface gravity data. 

The fast Fourier transformation has been applied to the computations. 


1. Introduction 

The formulas for the determination of the coefficients of the spherical harmonics of the earth's 
gravitational potential require the free-air anomalies to be given on a sphere, at least at a simp e 
See, e.g., ellipsoid or sea level. Thus we have to continue the free-air gravity anomalies 

downward to a sphere or an ellipsoid. 

The validity of the analytically downward continuation of the free-air gravity anomalies inside 
the earth is guaranteed by Runge’s theorem (Moritz, 1980, p. 67). Of course these gravity 
anomalies are not the original gravity anomalies inside the earth. The downward continuation 
gives a fictitious gravity anomaly on the ellipsoid that generates a disturbing potential T on and 

outside the earth. 

Moritz (1980) suggested that the free-air anomalies be continued to the point level. We can 
take also the ellipsoid as the reference surface and use this method to continue the gravi y 
anomalies down to the ellipsoid. 

Pellinen (1966) studied the methods for the determination of the coefficients of the spherical 
harmonics of the earth’s gravitational potential and he added a term, which can be easily 
transformed into terrain correction, to the free-air anomalies. 

In this paper we carried out some numerical investigation with the above mentioned methods. 
The terrain correction and the gradient solution were computed on a global basis. 


2. Mathematical formulation 

We continue the free-air anomalies downward to the ellipsoid by using the "gradient solution 
(Moritz, 1966, p. 68): 
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Ag = Ag + g , = Ag - h 


3Ag 

3h 


( 1 ) 


where Ag*, Ag are the gravity anomalies on the ellipsoid and on the earth's surface, respectively. 
The gradient of the gravity anomaly Ag is given by 


dAg R 


-hff 


Ag-Ag, 


dh 2 tc J H 3 q 


da 


( 2 ) 

where jlo = 2R sin vg/2, \jr is the angular distance between the current point and the computation 
point p, a is the unit sphere, and R is the mean radius of the earth. 


Pellinen (1966, p. 70) suggested that 
Ag* = Ag + G' 
where 


G =—{ f ( h ~ h p)( A g- A g P ) 

' JJ - Vo 


4 k “ ° 


da 


The relationship between the gravity anomaly Ag and the elevation h is assumed as 

Ag - a+ bh, b = 27tkp = 0. 1 1 mgal/meter 

a is a constant. 

Taking the plane approximation and using the assumption (5) we get 


gi — hpkp// -—dxdy 




(3) 
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( 6 ) 


G ~ C "I k P If TT^dy 


1 3 
0 


(7) 


w!th j^o = [(x-xp) 2 + (y-yp) 2 ] 1 / 2 , where x is the two-dimensional plane. 


109 


3. Computations of the Terrain Correction and the Gradient Solution on a global Basis 


The elevation data in 5' x 5' mean block values are available from the National Geophysical 
Data Center, Boulder, Colorado as ET0P05. This elevation data was used for the computation of 
the gi term and the terrain correction on a global basis. 

The integration region is taken as 1 5° in latitude extent and 30° in longitude. The boundary 
(or overlap) of the integration is taken 50 km. This satisfies the accuracy for most situations (Noe, 
1980). But the boundary is not large enough for the Himalaya Mountains. In the computations we 
took a 250 km boundary for the high mountain areas and a 50 km boundary for the flat areas. 

The gi term and the terrain correction are computed in 5' mean block values by using the FFT 
(fast Fourier transformation) technique. All computations are completed at the Instruction and 
Research Computer Center at Ohio State University. The CPU time required was about 10- 
second for each point on the IBM 3081 at OSU. 

The statistics of the gi term and the terrain correction are exhibited in Tables 1 and 2. 


Table 1. Statistics of the gi Term in 5' and 30' Mean Block Values 
Unit: mgal 


Block Size 

Mean Value 

Standard Dev. 

Max. Value 

Min. Value 

5' 

30' 

0.27 

0.27 

±2.56 

±1.54 

442.14 

45.08 

-78.88 

-10.47 

1° 

0.27 

±1.24 

25.52 

- 5.16 


Table 2. Statistics of the Terrain Correction in 5' and 30' Mean Block Values 
Unit: mgal 


Block Size 

Mean Value 

Standard Dev. 

Max. Value 

Min. Value 

5' 

0.23 

±1.01 

183.57 

0.6 

30' 

0.23 

±0.82 

25.24 

0.0 

1° 

0.23" 

±0.74 

17.77 

0.0 


It shows that the mean values of the gi term and of the terrain correction are almost the same. 
Of course, the gi term is larger and rougher (larger standard deviation, larger maximum values). 

After the computations of the gi term and terrain correction in 5' mean block values on a 

global basis, gi and the terrain correction were expanded in the spherical hamionics up to 180 th 
order. The RMS (root mean square) values of the degree variances of the gi term and the terrain 
correction are about 2 percent of the RMS values of the degree variance of the OSU86E gravity 
model. 

The total contributions of the gi term and the terrain correction on the geoid and the deflections 
of the vertical are shown in Table 4. 
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Table 4. RMS values of the corrections of the geoid undulation and the deflections of the 
vertical due to the gi and TC 


correction 

RMS ofSN 

RMS of 50 

gl 

0.7071 m 

0.1 129 secs 

TC 

0.7065 m 

0.0880 secs 


4. Conclusion 


The gi term and the terrain correction were computed in 5' mean block values on a global 
basis. The maximum of the gi term is 442 mgal located in the Himalaya Mountains. The 
maximum of the terrain correction is 184 mgal. The influence of the gi term and the terrain 
correction on the geoid undulation has been considered. It takes the order of 1 meter (RMS of 
correction of the geoid undulation). For the deflections of the vertical the RMS of the corrections 
of the gj term and the terrain correction is on the order 0.1". 
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